// Estimation du nombre pi par méthode de Monte Carlo

// clear;  // efface toutes les variables en mémoire.

// Scilab limite la taille possible des tableaux.
N = 1000000;

// Test s'il faut augmenter la taille de la pile mémoire des variables
st = stacksize();
if (8*N + 100000 > st(1)) then
  stacksize(8*N+100000);   
end

x = rand(N,1);  // = liste de  N  nombres aléatoires entre 0 et 1
y = rand(N,1);  // = liste de  N  nombres aléatoires entre 0 et 1

// Compte le nombre de couples de nombres (x(j), y(j)) tel que
// x(j)^2 + y(j)^2 < 1,  c'est-à-dire tel que le point
// de coordonnée  (x(j), y(j)) soit dans le cercle unité.
DansCercle = sum(x.^2 + y.^2 < 1);

disp("Approximation de Pi = " + string(4*DansCercle/N));
